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ABSTRACT 

We present an introduction to modern theories of interfacial fluctuations and the 
associated interfacial parameters: surface tension and surface stiflhess, as well as their 
interpretation within the capillary wave model. Transfer matrix spectrum properties due 
to fluctuation of an interface in a long-cylinder geometry are reviewed. The roughening 
transition and properties of rigid interfaces below the roughening temperature in 3fi 
lattice models are surveyed with emphasis on difl'erences in fluctuations and transfer 
matrix spectral properties of rigid vs. rough interfaces. 
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1. Interfaces at Phase Coexistence 

This article presents a review of various concepts and ideas in theories of inter- 
faces and their fluctuations. The presentation aims at providing a tutorial intro- 
duction emphasizing notation, nomenclature and uniformity of exposition of results 
accumulated in the field of Statistical Physics, on static interfacial fluctuation prop- 
erties. Thus, no completeness of material is attempted, nor is the list of literature 
complete or even balanced. Instead, general citation of review articles for further 
reference is provided (Refs. 1-12). 

Study of interfaces is an old topic dating back to first advances in Condensed 
Matter Physics. However, it grew immensely in scope and subject coverage re- 
cently, due to modern theoretical threatments with focus on fluctuation and phase- 
transition phenomena. Topics such as capillary waves, finite-size effects, wetting 
and surface phase transitions, roughening, connections with Polymer Physics and 
studies of membranes, etc., have each become a field of active research. The scope 
of the present review is limited to an overall introduction to selected topics, as well 
as summary of some recent results on the relations between interfacial fluctuations 
and transfer matrix spectral properties. As part of the conference Proceedings, 
and in fact, being a summary of a review talk presented, this survey avoids closely 
related topics covered in other presentations, notably, some finite-size effects. In- 
stead, emphasis was put on issues of relevance in numerical studies of interfacial 
properties. 
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Fig. 1. Interface in a narrow capillary and bounding a liquid drop. 

Static interfaces, or sharp domain walls, separate phases at coexistence. One 
has to assume that no gradual, continuous variation from one phase to another 
is possible, via a diffuse domain wall such as found, for instance, in 0(n) vector 
ferromagnets with continuous symmetry, n = 2, 3, . . ., and the associated soft-mode 
spin-wave excitations. Thus, interfaces separate distinct phases such as ± magne- 
tized (n = 1) Ising-model phases. In fact, walls and surfaces are all around us, and 
in the "zeroth-order approximation" that's all we actually see! However, thermally 
equilibrated (on experimental time scales) interfaces of the type Statistical Physics 
treats are less abundant. They can be observed between fluid phases, such as the 
gas-liquid meniscus in a narrow capillary, or surface of liquid drop on a wall; see 
Figure 1. The longest wave-length fluctuations are, however, typically cut off by 
gravity. True equilibrium crystal shapes can be also experimentally realized in some 
instances. 

Theoretically, the issue of how to define interfaces microscopically, and how to 
measure their properties, has encountered difficulties related to the fact that inter- 
facial fluctuation are rough, soft-mode, in many cases; details will be given later. 
Rough interface can be viewed as a sheet-like object, structureless and fluctuating, 
when it is "resolved" on a length scale large compared to the microscopic correla- 
tion lengths of the coexisting phases. This is illustrated by the broken-line box in 
Figure 2. A more microscopic resolution would reveal a region structured differ- 
ently from each of the two phases in contact, and forming a boundary between the 
phases; see the solid-line box in Figure 2. Theoretical description of the structure 
in this boundary region in terms of an "intrinsic" interfacial profile, e.g., density 
profile for a liquid-gas interface, is intimately connected with such box-resolution, 
so-called coarse-graining notions, on various length scales. 
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Fig. 2. Coarse-graining to reveal large-scale interfacial structure (broken- line box) 
and "intrinsic" microscopic structure (solid- line box), shown schematically. 

The concept of "intrinsic" vs. "large-scale" interfacial structure has never been 
well defined outside the framework of a given theory. There is no clear "snap- 
shot" definition of an interface. By examining a fixed system configuration, such as 
counting H — Ising bonds, identifying ± clusters, etc., one cannot define an inter- 
face unambiguously. There will always be overhangs, islands, clusters, which may 
be assigned to the phase-separating interface region, or equally well attributed to 
the single-phase fluctuation structure surrounding the interface. 

Thus, in order to define interface location and measure its properties, one appeals 
to external, essentially boundary-condition-based methods of fixing the interface. 
For instance, the interface in Figure 2 was secured by the ± boundary conditions at 
the upper and lower faces of the finite-size sample. If the sides are subject to periodic 
or free boundary conditions, then this interface will "float," oriented horizontally 
on the average. However, one can also pin the interfacial circumference by more 
restrictive boundary conditions. In Figure 3, the interface is forced to have average 
inclination angle 9 with respect to the horizontal axis of the sample, aligned with 
the axes of the underlying lattice (drawn schematically in the lower right corner). 
There are instances of interfaces generated spontaneously due to energy gain (such 
as in microemulsions) or entropy gain (for instance, in long-cylinder geometries 
with periodic or free boundary conditions in the transverse directions). However, 
we restrict our attention to interfaces forced in by boundary conditions. 

In a typical mixed-boundary-condition geometry there is an excess free-energy 
associated with the interface. For instance, in the geometries shown in Figures 2 
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Fig. 3. Interface with ends pinned by boundary conditions, in a 2id sample. Or cross-section of 
pinned-circumference interface of average area ^ j COS ^ in a f^-dimensional sample. 

and 3, the total free energy will be larger than in similar reference systems obtained 
by replacing all the — boundary conditions by +. The total free-energy increase 
will be proportional to the interfacial area. Let S denote the excess free energy per 
unit interfacial area, while F refer to the total free energy of the finite-size sample. 
We define 



L'^-'^l cose ' 



(1.1) 



where L'^''^ / cos 9 is the interfacial area, with ^ = for the case of Figure 2. The 
subscripts of the free energies indicate the boundary conditions, as just described. 

The interfacial free energy density, 11{T,9, L, . . .), depends on system sizes, on 
shape ratios and boundary conditions. Most of these dependencies are, however, 
corrections typically of order 1/L or smaller. The leading, order 1, contribution in 
the thermodynamic limit of infinite sample, gives the macroscopic surface tension. 
Finite-size corrections in some geometries are surveyed in other contributions to 
these Proceedings. 

It is convenient to introduce the reduced quantity, of units (area)~^. 



a = T,/kT 



(1.2) 
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An important observation is that interfaces usually separate phases which are non- 
critical and therefore microscopic fluctuations are not necessarily isotropic. Thus, 
not only the temperature dependence but also the angle dependence of the surface 
tension survives in the thermodynamic limit. The angular dependence is with re- 
spect to the underlying lattice structure, so that it is not present for fluids but plays 
an important role when at least one of the coexisting phases is solid. For model 
lattice systems used in numerical simulations, the finite-size samples are usually 
oriented along a symmetric lattice axis, as marked schematically in Figure 3. In d 
dimensions, (d — 1) angular variables are needed to specify the angular dependence 
of cr. 

2. Surface Tension and Surface Stiffness Coefficient 

In this section we survey properties of the reduced bulk (i.e., infinite-system 
limiting) surface tension, (t{T,9). We use the 2d notation with one angle, and 
whenever appropriate, quote results in the form relevant for rough interfaces. Both 
the concept of rough interfaces and the transition to the rigid-interface behavior at 
low temperatures in lattice 3d models will be further explained in later sections. 

Interfaces are termed "rough" when due to gain in entropy they are not pinned 
by the underlying lattice structure. Their fluctuations are long-ranged, soft-mode. 
They can therefore pass locally at any angle with respect to the lattice axes. The 
angular dependence is smooth, analytic. 



As mentioned, in numerical studies the angle is always measured with respect to a 
symmetric axis, so that the term linear in 9 vanishes. 

In the geometry of Figure 3, let us consider the interfacial free energy per unit 
projected, horizontal area, L^~^ , cf. (1.1), 



where we defined an important quantity 



termed the surface stiffness coefficient. This nomenclature will become clear in the 
next section, when we describe the capillary wave model. 

In two dimensions, interfaces are always rough, for < T < Tc, where we use 
the Ising-model notation for the ± phase coexistence in zero applied field below 
the Curie temperature Tc. As shown in Figure 4, the 2d surface tension vanishes 
according to 




(2.1) 




(2.2) 



(7 ~ (Te - T) 



(2.4) 
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Fig. 4. Surface tension of — 2 systems. 



divergence 



kTK 



rigid 




Fig. 5. Surface stiifness coefficient of — 2 systems. 



as T —f T~ , with the exponent values typically /i = 1 or smaller (infinite slope), 
depending on a system universality class. At T = 0, the interfaces are rigid even in 
2d. This concept will be described in detail in Section 5 below. For now it suffices 
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Fig. 6. Surface tension of = 3 systems. 

to note that rigid interfaces do follow lattice-axis directions microscopically, pinned 
by energy constraints at low temperatures when entropy gain plays no important 
role. The energy cost of the interface is finite at low temperatures; the unreduced 
quantity kTcr = S, shown in Figure 4, approaches a finite value at T = 0. 

The surface stiffness coefficient in 2d is shown schematically in Figure 5. Near 
Tc, where the fluctuations become isotropic, the added second-derivative term in 
(2.3) is asymptotically smaller than the first term. The surface stiffness simply 
follows the surface tension. However, as T ^ 0, the second term dominates. In 
the expression fcT/t, plotted in Figure 5, the surface-tension contribution is finite 
as T ^ 0. However, the second-angular-derivative contribution actually diverges 
according to ex.p (const /kT), as T —f 0, where the constant is system-dependent. 

The surface tension of the 3d systems, see Figure 6, behaves similarly to the 
2d case. Near Tc, (2.4) applies, with typical ji values in the range 1.0+ to 1.5, 
depending on the universality class of the critical point at Tc, at which the phases 
loose their distinction and the interfacial structure merges with the rest of the 
homogeneous system. However, 3d interfaces of lattice (crystalline) systems only, 
become rigid for < T < Tr, where the roughening temperature is of order Tc/2. 
The surface tension has an essential singularity at T = Tr. 

The surface stiffness coefficient in 3d is shown schematically in Figure 7. As 
will be further described later, in Sections 3-5, it is defined only for Tr < T < Tc- 
Indeed, it turns out that in the "rigid" regime we have the angular dependence in 
the form cr(T, \9\), where 9 is measured with respect to one of the lattice axis, so 
that (2.1)-(2.3) cannot be used. For a more detailed description of the behavior 
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Fig. 7. Surface stiffness coefficient of lattice d — 3 systems. 

near Tr the reader is addressed to a review on the roughening transition (Ref. 9) 
which is believed to belong to the Kosterlitz-Thouless universality class of d = 2, 
0(2)-symmetric bulk XY systems. 

3. Capillary Wave Model 

In this section we describe the capillary wave model of the longest wave-length, 
"soft-mode" fluctuations of rough interfaces. Our formulation will be in the 2d 
notation, while extensions to d > 2 will be described as appropriate. On large 
length scales, rough interfaces look sheet-like and can be approximately described 
by a single valued height function h(x, . . .), where the unspecified arguments refer 
to dependence on other transverse coordinates, say y, z, . . ., for d > 2. 

In the capillary wave model, the energy of the fluctuating interface is approxi- 
mated by the "macroscopic" energy obtained by intergrating the angle-dependent 
surface tension over the surface defined by the height function h. The 2d definitions 
are illustrated in Figure 8. We note that 

fdh\ , , 

Q = arctan I ^ I j (^-1) 

and so the energy is given by 



E = kT J dx (T[arctan (dh/dx)] \J 1 + {dh/dxf 
The square root in (3.2) comes from the element of length. 



(3.2) 
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Fig. 8. Capillary wave interface description. 



d£ = ^{dhf + {dxf . (3.3) 

Thus, for d > 2 the derivative term under the square root will be replaced by 
(dh/dx)^ + [dhj dyY + . . ., with (d — 1) transverse-coordinate derivatives. The 
angular dependence in d > 2 will also be more complicated, as already mentioned. 

The single-sheet description ignoring overhangs, bubbles, and generally assum- 
ing that all the "intrinsic" interface structure has been absorbed in the macroscopic 
surface tension values, is of course approximate. In fact, to make the resulting 
model tractable, one further assumes that the interface remains largely flat. Thus, 
the derivatives are small. Specifically, the integrand in (3.2) can be expanded in 
(dh/dx), keeping the second-order terms. 

The prime here denotes the angular derivative. Note the emergence of the surface 
stiffness coefficient as the prefactor of the quadratic term, — the origin of the term 
"surface stiffness." 

Integration of the first term yields the constant reference energy proportional 
to the interface length. The second term only contributes at the end points, the 
precise specification of which is of no interest. Only the quadratic term depends on 
the local fluctuating variable h(x). The resulting expression 
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CW 




(3.5) 



is the standard capillary wave energy, Ecw ■, in 2d. 

In higher dimensions, for general interface orientations, expansion to the second 
order in derivatives with respect to the (d — 1) transverse coordinates will produce a 
mixed-derivative (second order) form due to angular dependence. In fact, a rotation 
to principal-axes coordinates, (a;, y, z, . . .) (X, . . .), is required to obtain a 
diagonal form 



However, in most numerical applications the interfacial orientation is symmetric 
with respect to the underlying lattice structure. Thus, the original coordinates can 
be used, and all the stiffness coefficients are equal. 

Strictly speaking, the capillary wave energies (Hamiltonians) (3.5), (3.6), must 
be supplemented by cutoff prescriptions to have finite expressions for thermody- 
namic average quantities, when these energies are used in Boltzmann factors. Oth- 
erwise, these Gaussian models are exactly solvable. Specifically, in 2d the capillary 
fluctuations are indeed strong. Calculations with the Hamiltonian (3.5) suggest 
that a 2d interface "viewed" on the length scale L (for instance, pinned by its ends 
at separation L), fluctuates over vertical (/i-direction, see Figure 8) span of order 
\/L. Detailed Gaussian calculations for all d, and references to literature, can be 
found for instance in Refs. 13, 3. 



The case d = 3 turns out to be marginal: the interface only spans ~ V In L in the 

vertical direction, when viewed on length scales L (in both transverse directions). 
For d > 3, the capillary span is of order 1. Thus, the model must be applied with 
great caution. When capillary fluctuations are comparable to other fluctuations 
in the system, there is no obvious reason to single them out. Some indirect evi- 
dence for characteristically Gaussian, capillary fluctuations may come from studies 
of certain finite-size effects. However, the form of these has been a subject of a 
recent controversy which will not be reviewed here (but see other articles in these 
Proceedings). 

4. Transfer Matrix Method and Capillary Wave Fluctuations 

Interfacial parameters, cr, /t, etc., can be measured by various methods. As 
already mentioned the most direct approach is based on comparing a system with 
an interface with another system where the boundary conditions do not impose any 
interfacial structure. While finite-size effects are usually corrections to the leading 
"bulk" results, they can also be used to measure interfacial parameters by numerical 
techniques. 




(3.6) 
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In this section we describe yet another method of measuring the interfacial 
properties, based of spectral features of the transfer matrix (Refs. 14-15). No review 
of the transfer matrix method will be presented here (see Ref. 16-18). Instead, we 
offer the following general outline. The transfer matrix method is based on summing 
up over the system configurations by adding consecutive slices of size 6 x ^4, in the 
cylindrical geometry of (d— l)-dimensional cross-section of area A, and axial length 
L|| . Thus, the partition function is calculated by summing up over configurations 
of L^^/b slices. In numerical calculations b is usually just one lattice spacing, while 
the cross-section shape is square or other lattice-aligned shape. 

Boltzmann factors incorporating the interaction energy of the added slice with 
the previous one, as well as its internal energy (to have a symmetric transfer matrix 
the latter is usually replaced by the mean over two slices), are arranged in a matrix 
labeled by slice configurations. Calculation of the partition function by summing 
over system configurations then reduces to matrix multiplication, with summation 
over each matrix index corresponding to summing over the configurations of one 
slice of axial length b. 

Of particular interest are the eigenvalues of this "transfer matrix," 



Ao > Ai > A2 > . . . , (4.1) 

some of which can be multiply-degenerate (but not Aq). The reduced (per kT) 
free-energy density of the infinite system, L» = 00, is given by 



/=-(lnAo)/(6^). (4.2) 

This free energy of course depends on the cross-section size, shape, boundary con- 
ditions, etc., as corrections. It approaches the bulk free-energy density, /(oo), as 
^ ^ 00. 

Correlation lengths (proportional to the inverse energy gaps in the Quantum- 
Hamiltonian nomenclature popular in Particle Physics) in the long-cylinder geome- 
try, L|| = 00, describe the exponential decay of correlations at large distances along 
the cylinder. They are given by 

= ln(A!/A,) ' 

with i = 1,2,.... 

Correlation lengths associated with the longest-length-scale fluctuations will be 
the largest, i.e., they will have the lowest indices j in (4.3). Generally, the number 
of slice configurations increases exponentially with the cross-section size. Numeri- 
cal diagonalization of transfer matrices is thus a large-scale computational problem. 
However, theoretically one can in some cases express the largest correlation lengths 
explicitly in terms of few parameters. This description relies on that fluctuations 
in the long-cylinder geometry are effectively one-dimensional on large length scales. 
Thus, an effective, few-parameter parametrization of the transfer matrix can be 
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found (Refs. 19, 15) which reproduces correctly the asymptotic form of the largest 
correlation lengths. As a result, parameters determining large-scale fluctuations, 
which for systems phase-separated in the bulk, usually turn out to be various in- 
terfacial properties, can be probed directly by transfer-matrix calculations. We 
describe one such application in this section. 

Let us consider the long-strip 2d geometry shown in Figure 8, with ± boundary 
conditions ensuring that an interface is running along the strip. The height of 
the strip in the transverse direction will be denoted by H, and we assume that 
it spans the distance from h = to h = H . (The axes in Figure 8 were shifted 
for clarity.) In the effective, phenomenological description focusing only on the 
large-scale interfacial fluctuations which are rough (we assume < T < Tc), the 
configuration is specified by the height variable h. The effective transfer matrix 
becomes an operator "labeled" by h, in the following sense. The vector space on 
which it acts is that of functions tpih), and the boundary conditions which essentially 
mean that the interface cannot cross the walls, can be shown to correspond to 



V'(O) = V'(-ff) = 



(4.4) 



The effective transfer operator acts on these functions, and it turns out that the 
/i(a;)-related fluctuation part of the interfacial energy, singled out in relation (3.5), 
is represented in this language by the operator 



b d? 



(4.5) 



We will not describe in detail the derivation of this result but instead emphasize its 
physical meaning and implications. The "free-field Hamiltonian" appearing in the 
exponent of (4.5) naturally corresponds to the fact that interfacial fluctuations are 
diffusion-like in the transverse direction, within the capillary wave model, with the 
diffusion rate inversely proportional to the surface stiffness coefficient k. 

As in the derivation of (3.5), here there are also contributions which do not arise 
from the fluctuations included in the capillary wave description. These account for 
the overall free energy density of the bulk, fcT/(oo), of the interface (flat), kTa 
(measured per unit length), as well as two ±-wall free-energy terms, equal due to 
the ± symmetry. The latter will be denoted by kTw per unit axial length. Thus, 
the effective transfer operator takes the form 



exp 



-bHf (oo) -b{(T + 2w) 



b 

2k dh^ 



(4.6) 



The eigenvalues are all nondegenerate. 



A, = exp 



-bHf (oo) - b{a+2w) - b 



^^(i + 1)^ 
2kH^ 



(4.7) 



The free-energy relation (4.2) with A = H then gives 
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f{H) = /(oo) + ^ + ^ + « {H-') . (4.8) 

Thus, we derived the leading-order capillary correction to the free energy in the 
strip geometry. While the first two terms are the standard reference bulk, surface, 
and interface free energies, the fluctuation contribution of order is determined 
by the surface stiffness coefficient. Note that the higher-order terms are not known 
generally (the lower-case o denotes "less than order of"). They need not be fully 
capillary wave dominated. 

The correlation lengths follow from (4.3) as 

This form of course applies only to the leading correlation lengths, j = 1, 2, . . .; the 
precise criterion is described elsewhere (Ref. 15). The largest correlation lengths 
thus provide a direct measure of k. 

Various other results and observations derivable from the effective transfer- 
matrix formulation are not reviewed here. Let us only mention that the "ground 
state" wave function tpoih) is real, non-negative, and provides the probability dis- 
tribution of finding the interface at h. There are other differences between the 
quantum mechanical time evolution formulation, in "Euclidian time," where the 
wave functions yield probability amplitudes, and the transfer matrix formulation 
in Statistical Mechanics. Thus, caution and care must be exercised when methods 
and results are taken over from one field to another. 

In numerical studies, systems with antiperiodic boundary conditions are fre- 
quently preferred to those with walls because the former preserve some of the 
translational-invariance features, and there are no wall energies, it; = 0. The trans- 
fer matrix formulation can be easily extended to antiperiodic boundary conditions. 
However, it turns out that these boundary conditions effectively induce a topology 
of mod(2ff) in the strip: interfacial height h can be reached with + above, — below 
the interface, or with — above and + below. Two rotations of the interface around 
the system, by vertical displacement, are required to repeat the same vertical-slice 
configuration of the ± phases. 

For antiperiodic boundary conditions, the relation (4.4) is replaced by 

■>i){h) = ■>i){h + 2H) . (4.10) 
The results (4.8) and (4.9) are then replaced by 

f{H) = /(oo) + |: + O [exp {-H/Uiu)] , (4.11) 
^.>o=^ + «(^^) , (4.12) 
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H 




Fig. 9. Interface in a long-cylinder geometry in 3d. 

where the correction in (4.11) is believed to be exponentially small. Other differences 
include double-degeneracy of Aj>0) etc. 

Exact calculations for 2d Ising-model strips yield expressions for the correlation 
lengths (Refs. 20-23) which check the capillary wave model predictions and also 
verify the limits of validity of the capillary approximation which will be briefly 
discussed at the end of this section. 

We now turn to the 3d case, in the rough-interface regime, Tr < T < Tc- The 
± long-cylinder geometry is shown schematically in Figure 9. The system size in 
the direction perpendicular to the (h, x) plane (marked by y in the figure) will 
be denoted by P. We assume that the boundary conditions along the third, y- 
axis direction are periodic or free, so that the only interface-imposing boundary 
conditions remain at the top and bottom walls (at h = 0, H), as in the 2d case. We 
also assume that P/H = 0(1). 

The longest wave-length capillary modes then "propagate" along the cylinder. 
They correspond to fluctuations that leave the interface relatively flat along the y 
direction; recall that in regions of axial size Aa; = 0(P) the interface is only spread 
O (y^og P) . However, the interface can fluctuate up and down as a whole, on large 
s-direction length scales which are precisely the largest correlation lengths that we 
want to estimate. However, these fluctuations are effectively two-dimensional, with 
the effective stiffness coefficient Pkx- 

We adopt the notation suggested by (3.6): the 3d stiffness coefficient will be 
denoted by Kx- Then the effective value for the 2d-like fluctuations IS K' — P tx/^ . 
and we can use the previously derived leading-order expressions (4.8), (4.9), (4.11), 
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(4.12) with this identification. Note that the units of k are (length) -(<*-i) so that 
the effective 2d stiffness coefficient has the correct units. 

Recent numerical study (Ref. 24) used the capillary wave model results described 
here to estimate the surface stiffness coefficient of the 3d Ising model in the rough- 
interface regime. Both the exact 2d Ising results mentioned earlier and the numerical 
3d studies indicate that on approach to the "dangerous" points in the phase diagram, 
T = Tc and, in 3d, T = Tr, the expressions derived above must be used with 
caution. Specifically, for T —f T~ , the forms proposed will "cross-over" to the 
full finite-size scaling description. The corrections in inverse powers of H in fact 
become functions of ^luik/H, and the capillary wave expressions break down as 
^buik diverges at criticality. 

At the roughening temperature, the finite-size scaling forms are not well un- 
derstood, and in fact it is not even clear if they exist. Thus, we have no theory 
for the form of cross-over from the rough-interface capillary wave regime to the 
rigid-interface regime (see the next section), in finite-size systems. 

5. Rigid Three-Dimensional Interfaces at Low Temperatures 

In order to complete our discussion of interfacial fluctuations, we now address 
fluctuation properties of interfaces in the low-temperature regime of 3d lattice mod- 
els (or crystalline-phase systems) which, as already pointed out are "rigid." More 
precisely, at low temperatures the interfaces prefer to align with the lattice struc- 
ture due to energy gain. Thus, an interface which on the average is oriented along 
one of the underlying lattice axes, will follow this direction not just on the average 
but also microscopically. Of course, for any nonzero temperature the interface will 
be "dressed" by thermal fluctuation features such as pits, overhangs, bends, etc. 
However, these are no soft-mode excitations as for rough interfaces. Rather, these 
thermal fluctuations are on length scales comparable to the bulk correlation lengths 
of the coexisting phases. 

The two interfacial fluctuation regimes are separated by a sharp phase transition, 
the roughening transition (Refs. 9), which is in the Kosterlitz-Thouless universality 
class. However, our discussion will be focused on fluctuations away from Tr. In 
fact, we assume that the temperature is fixed in < T < Tr, and the system 
sizes are large enough so the the "critical-like" aspects of the roughening as a phase 
transition have no effect. (See also note at the end of the preceding section.) 

The nature of the rigid-interface fluctuations is best revealed by considering the 
so-called "vicinal" interfaces which are inclined a small angle 9 with respect to the 
preferred symmetric lattice direction. Such as interface is schematically shown in 
the right panel of Figure 10. As mentioned, rigid interfaces tend to follow lattice 
directions. For simplicity let us assume that the lattice is cubic, of spacing o, as 
shown schematically in the lower right corner of the picture. The interface shown 
in Figure 10 is forced by boundary conditions to raise (Ah/a) lattice spacings when 
it spans the sample length L. Rough interfaces follow their average inclination, 
see the left panel in Figure 10, accompanied by large, capillary fluctuations on 
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Fig. 10. Rough vs. rigid interfaces in 3d. 

d = 3 




Fig. 11. Step free energy of rigid interfaces in 3d. 

scales exceeding the lattice spacing o. However, rigid interfaces actually achieve the 
appropriate inclination angle by having (Ah/ a) microscopic steps. 

Each step has excess free energy proportional to its length, M , where we assume 
that the sample dimensions in the horizontal cross-section are Lx M , see Figure 10. 
Let s{T) denote this excess free-energy per unit length (of a step) and per kT. Ac- 
tually, s{T) is defined only with respect to symmetric lattice directions, and in fact 
not only its T-dependence may vary for different directions, but the roughening tem- 
perature, Tr, itself may differ for various lattice orientations as well. The behavior 
of this step free energy is shown schematically in Figure 11. At low temperatures. 
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kTs is a finite quantity. As T —f , the step free energy vanishes with an es- 
sential singularity. Near the roughening transition the length scale diverges 
and plays the role of a correlation length similar to the correlation length at the 
Kosterlitz-Thouless transition in 2d XY systems on the high-temperature side of 
that transition (Refs. 9, 25-26). For T > Tr, formally s = 0. 

The total free-energy cost of the rigid interface shown in Figure 10 is given by 

I\F = kTa{0)LM + kTsM{\Ah\/a) + ... . (5.1) 
Since the inclination angle 9 is 

e = arctan (AV-Z^) , (5.2) 
the angular-dependent surface tension takes the form 

Note for comparison that the leading angle-dependent free energy term for rough 
interfaces in symmetric lattice directions is quadratic in 9, as suggested by (2.1) 
with c'(0) = 0. However, in the rigid-interface regime the ^-dependence is no 
longer analytic, and the leading, order |^| free-energy correction is determined by 
the step free energy and is explicitly lattice-spacing dependent. 

Of cause the above arguments only apply for very small angles, when the steps 
are well-separated. The correction terms shown by ... in (5.1) and (5.3) will be 
contributed by the interplay of two effects. Firstly, at any nonzero temperature 
the steps, as well as the flat interface portions, see Figure 10, will be "dressed" by 
thermal-fluctuation features. Those on scales of the correlation lengths of the two 
coexisting phases can be considered absorbed in the definition of s{T). However, 
the steps when viewed "from above" in Figure 10, will not be just straight lines of 
length M. They can have kink-like shifts along the L-direction and will therefore 
behave random-walk like. Even though "kinks" cost energy, and as a result the 
effective two-dimensional stiffness of these soft-mode-like step fluctuations will be 
large, it will definitely not be infinite. 

Secondly, for finite angles 9, the step wandering just described will contribute 
a long-range term (Refs. 9) in step-step interactions. Short-range step-step in- 
teractions also arise from energy and entropy changes associated with step con- 
tacts. These interactions determine the corrections in (5.1), (5.3). Step interac- 
tion and fluctuation phenomena are accounted for semiquantitatively within the 
teiiace-ledge-kink model of rigid-interface fluctuations. The reader should consult 
the more specialized review (Ref. 9) for details, further literature citations, and for 
a description of how the large-angle variation of the interfacial free energy actually 
determines crystal shapes. 

Direct estimates of the step free energy have been reported in Monte Carlo 
studies of 3d Ising models with Ah = o (Refs. 25-26). In the rest of this section we 
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Fig. 12. Directed-walk model of rigid interface fluctuations in 3d cylinders. 

describe how the transfer matrix spectrum is changed in the rigid-interface regime. 
We find that the spectral gaps (inverse correlation lengths) are determined by s{T) 
and thus provide a direct measure of this interfacial property. 

In the long-cylinder geometry, of cross-section size H x P, see Figure 9, we 
now consider the regime of < T < Tr. The interface will no longer be rough 
as in Figure 9. Instead, its will be microscopically flat over large length scales, 
measured by the leading correlation lengths ^j. The dominant fluctuations will 
be well-separated steps (we will see later that ^ H,P). Such fluctuations can 
be modeled by a solid-on-solid model, of steps (forward, up and down) forming a 
directed random walk, as illustrated in Figure 12. Such models have been used 
extensively for two-dimensional interfaces (Refs. 3, 5, 11). 

However, as a model of the "cross-sectional" form of the rigid 3d interface the 
directed walk has a simplifying feature. Indeed, the transverse, up and down, steps 
of length o are extremely improbable. Their free-energy cost is represented by the 
exponentially small (in the system size P) Boltzmann factor 



exp(-Ps + ...) . (5.4) 

The Boltzmann factor per each horizontal step of length b, — the slice size, — is 
given by the expression that incorporates the nonfluctuation, "reference" free-energy 
contributions. 



exp [-bH Pf (oo) - bPa + . . .] . (5.5) 
The correction terms in (5.4) and (5.5), marked depend on finite-size effects 
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associated with the boundary conditions in the H and P directions; they will not 
be described here (see Ref. 15). 

Note that the weight of the transverse steps is small as compared to the "no 
step" Boltzmann factor, 1. Longitudinal weights are also numerically small, but for 
directed walks their number (per unit system length) is not a fluctuating quantity. 
Due to small weights (5.4), one can show that the results do not depend on the 
details of the walk model chosen. For instance, they are the same if only single or 
any number of transverse steps are allowed at each fixed s-coordinate value (these 
are known as, respectively, restricted and unrestricted 2d solid-on-solid models). 

The effective transfer matrix acts now in the space determined by the allowed 
/i-coordinate values which are no longer continuous as for rough interfaces. Indeed, 
as illustrated in Figure 12, the allowed h values will be discretized in steps of o. 
Thus, the interface can have O (H/a) height values. The effective transfer matrix 
acts on vectors with this number of elements and therefore has O (H/a) eigenvalues. 
Explicit calculations are rather technical, and we only summarize here the resulting 
expressions for the correlation lengths. 

The eigenvalues of the effective transfer matrix form a distinct multiplet asso- 
ciated with the rigid-interface fluctuations. Thus, the leading O (H/a) eigenvalues 
are in fact split in a very narrow band, ~ exp(— Ps), and they are well separated (by 
spectral gaps of order 1) from the rest of the spectrum (of the full transfer matrix). 
For ± boundary conditions, shown in Figure 12, we find the expression, 

-1 

(5.6) 

where the eigenvalues are labeled by j = 0, 1, . . ., (H/a) — 1. Thus, only the (H/a) 
largest, nondegenerate eigenvalues are determined by s{T) via the leading, expo- 
nential size-dependence of the (inverse) gaps. 

The prefactor quantities, 7...(T) and the power w, are boundary condition and 
shape-dependent. The power w has been related to finite-size corrections due to 
the finite extent P. For, respectively, periodic and free boundary conditions in the 
P direction, the predicted values are w = ^ and 0. However, these results are 
essentially just conjectures. 

For antiperiodic (op) boundary conditions, the results are 

^,>o = 6T.,(T)P-e^»(^)(^sin|^y' , (5.7) 

where now the index covers the range j = 0, 1, . . ., H/a, and furthermore all but 
the j = 0, H/a eigenvalues are doubly-degenerate. The power w is the same as in 
(5.6). No detailed predictions are available for the functions 7...(T). 

Numerical studies of the 3d Ising model below Tr (Ref. 24) checked the leading- 
order step-free-energy dependence, ~ exp(Ps). The expressions surveyed here work 
well except on approach to the roughening temperature Tr. In this limit the sizes 
needed to observe the asymptotic regime of validity of the "microscopic" expressions 



^,>o = 6T±(T)P™e^^(^) 



. 7ro(i + 2) 
sin — ; sin ■ 



2(H + a) 2(H + a) 
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(5.6), (5.7), etc., become large and cannot be reached in numerical simulations. As 
emphasized in Section 4, there is presently no cross-over theory to describe the 
behavior near Tr. 
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